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Abstract 

The potential of populations to evolve in response to ongoing climate change is 
partly conditioned by the presence of heritable genetic variation in relevant 
physiological traits. Recent research suggests that Drosophila melanogaster exhib- 
its negligible heritability, hence little evolutionary potential in heat tolerance 
when measured under slow heating rates that presumably mimic conditions in 
nature. Here, we study the effects of directional selection for increased heat 
tolerance using Drosophila as a model system. We combine a physiological 
model to simulate thermal tolerance assays with multilocus models for quanti- 
tative traits. Our simulations show that, whereas the evolutionary response of 
the genetically determined upper thermal limit (CTmax) is independent of 
methodological context, the response in knockdown temperatures varies with 
measurement protocol and is substantially (up to 50%) lower than for CTmax. 
Realized heritabilities of knockdown temperature may grossly underestimate the 
true heritability of CTmax. For instance, assuming that the true heritability of 
CTmax in the base population is h 2 = 0.25, realized heritabilities of knockdown 
temperature are around 0.08-0.16 depending on heating rate. These effects are 
higher in slow heating assays, suggesting that flawed methodology might 
explain the apparently limited evolutionary potential of cosmopolitan D. mela- 
nogaster. 
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Introduction 

The ability to adapt and tolerate ongoing rising tempera- 
tures depends to a large extent on organismal plasticity 
and the evolutionary potential of populations (Helmuth 
et al. 2005; Ghalambor et al. 2007; Angilletta 2009). In 
terrestrial ectotherms, the limited evidence suggests that 
many tropical and desert species live at temperatures near 
their upper thermal limits (Deutsch et al. 2008; Huey 
et al. 2009), resulting in mounting pressure to understand 
when evolutionary responses may counter rapid climate 
change and how to quantitate this evolutionary potential 
(Parmesan 2006; Skelly et al. 2007; Visser 2008; Hoff- 
mann and Sgro 2011). Current information appears to 
indicate that upper thermal limits are weakly correlated 
with latitude in terrestrial ectotherms (Addo-Bediako 



et al. 2000; Sunday et al. 2011), and recent studies have 
questioned the ecological relevance of artificial selection 
experiments that have successfully increased heat toler- 
ance in Drosophila because flies were placed acutely at 
stressful temperatures or subjected to a fast heating rate 
(e.g., McColl et al. 1996; Gilchrist and Huey 1999), which 
may overestimate species tolerance limits (Chown et al. 
2010; Hoffmann 2010). At slow and presumably more 
realistic heating rates, heat tolerance is substantially lower 
and its additive genetic variance, and consequently nar- 
row-sense heritability (henceforth simply "heritability"), 
are almost negligible (Terblanche et al. 2007; Chown 
et al. 2009; Peck et al. 2009; Mitchell and Hoffmann 
2010). These findings lead Mitchell and Hoffmann (2010, 
p. 699) to conclude (our addition between brackets): If a 
highly adaptable species like D. melanogaster which exhibits 
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generally high heritability estimates for all quantitative 
traits has problems mounting evolutionary responses [for 
upper thermal limits], where does this leave other species 
whose adaptive potential might be curtailed due to small 
population size or a history of selection? 

Even though the presence of evolutionary limits is cer- 
tainly a matter of concern, it is debatable whether limited 
evolutionary potential ultimately accounts for many 
empirical results. It is becoming increasingly evident that 
methodology can have a greater impact on estimates of 
upper critical thermal limits (CT max , defined as "the max- 
imum temperature that an organism might potentially 
tolerate given its physiological condition in the absence of 
any other hazard;" Santos et al. 2011) than the evolution- 
ary or acclimatory responses that researchers aim to study 
(Lutterschmidt and Hutchinson 1997; Terblanche et al. 
2007; Chown et al. 2009; Rezende et al. 2011; Santos 
et al. 2011; Ribeiro et al. 2012). This raises the possibility 
that numerous reports on evolutionary limits might 
ultimately reflect measurement artifacts because the 
magnitude of a thermal challenge during an assay 
depends on both the temperature and the duration of 
exposure to this temperature (Hochachka and Somero 
2002, p. 331), and an organism's physiology as well as its 
probability to survive any given thermal challenge vary in 
time. Taking these issues into account, we have recently 



positively correlated with AT (Elliott et al. 1994; Mora 
and Maya 2006; Terblanche et al. 2007; Chown et al. 
2009; Peck et al. 2009). Rezende et al. (2011) also demon- 
strated that ramping rates affect the additive genetic and 
residual variances of CT max in opposite directions, with 
slow ramping rates decreasing the genetic variance, but 
increasing the residual variance. However, in their model, 
the residual variance arises from individual differences in 
metabolism, which consume energy and water resources 
at rates that are unrelated to genetic differences in CT max . 

The theoretical approach in Rezende et al. (2011) is a 
simplified account of what happens in heat resistance 
assays, and was recently expanded by incorporating a sur- 
vival probability function that varies with temperature 
(Santos et al. 2011). In this model, knockdown tempera- 
tures involve a deterministic component, as in Rezende 
et al. (2011), and a stochastic component that reflects the 
time-dependent cumulative probability of collapse as tem- 
perature approaches CT max (see Appendix SI). As CT max 
corresponds by definition to the upper physiological limit, 
knockdown temperatures will always be biased 
downwards with respect to this parameter, and this is the 
primary reason why empirical knockdown temperatures 
should not be equated with CT max . The key ingredients of 
Santos et al.'s (2011) model are encapsulated in the 
following equation (eqn 10 in their paper): 



X f" if Tb < Threshold 



^\ ^rT^Tr^ n- >< if r » * r "ta. 

[{[CT m „((,)-r b ]/[CT m „((,)+r b ]} J 



developed a theoretical framework that adequately repro- 
duces the impact of methodological protocol on empirical 
measurements of CT max in D. melanogaster and explains 
why measurements of CT max obtained with different 
methods are often uncorrelated or exhibit contrasting 
heritability estimates (Santos et al. 2011). Here, we study 
the effects of directional selection for increasing heat tol- 
erance with a theoretical approach that combines this 
framework with multilocus models for quantitative traits. 

Background 

Rezende et al. (2011) postulated that CT max can change 
due to acclimation and resource depletion (or fatigue) 
during the course of a ramping assay, in which tempera- 
ture increases from an initial temperature T 0 at a rate AT 
(°C /min) until individuals succumb to heat stress. The 
total amount of time under heat stress increases in slow 
ramping assays, which lowers CT max as resources are 
depleted and explains why heat tolerance is often 



where p(f;) is the probability of any given individual 
surviving to time f ; (the time interval is 1 min) as a func- 
tion of body temperature T b ; Budget(t 0 ) is its total 
reserves at t 0 , which is depleted during the course of the 
experiment at a rate determined by metabolism; EC(ti) is 
the amount of resources consumed at time f,-; Threshold is 
the temperature at which or above which the individual 
is under thermal stress; and k and a are constants. The 
model can accurately replicate survival times when flies 
are assayed for desiccation resistance (Tb < Threshold! 
Fig. 1 in Santos et al. 2011) or are subjected to different 
types of thermal stress (T b > T thresho i d ; Fig. 3 in Santos 
et al. 2011). 

Whether those individuals with the highest knockdown 
resistance correspond to those with the highest CT max will 
depend on the amount of noise introduced by stochasticity, 
which is expected to be greater in heat tolerance assays 
that use slow heating rates (Santos et al. 2011) that are 
presumably "ecologically realistic" (Chown et al. 2009; 
Mitchell and Hoffmann 2010; Overgaard et al. 2011a,b; 
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Figure 1. Sample numerical results from simulation model 1 assuming additive and equal allele effects for CT max . The census size at each 
generation was N = 5, 000 flies. They were subjected to 50 generations of up-selection for knockdown temperature (the top 20% of each sex 
was retained) using a fast ramping protocol with 7" 0 = 28°C and AT = 0.5°C/min. CT max was controlled by 20 diallelic loci on the same 
chromosome, with recombination frequency between adjacent loci r = 0.25 in females. Allelic frequencies in the base population ranged from 
p = 0.1 to p = 0.5 for alleles Y, which all have additive effects 0.1°C. Heritability of CT max was h 2 = 0.25 in the base population. Panel A plots 
the increase in knockdown temperature and its realized heritability, estimated by regressing the response to selection against the cumulated 
selection differential over the first 12 generations of selection. Panel B plots the increase in CT max , which was 20% higher relative to the increase 
in knockdown temperature. Panel C plots the frequency changes of alleles Y increasing CT max , which eventually reached fixation. Panel D plots 
the total genotypic variance V G together with its causal components. The genie variance component V g initially increased (approximately up to 
generation 8) as a consequence of the changes in allele frequencies, but the genotypic variance V c steadily decreased because of linkage 
disequilibrium (D L ). Plots for CT max , allele frequencies, and variance components are framed in shadow because their responses to selection for 
knockdown temperature are hidden to the experimentalist. For results of a similar model relaxing the assumptions of equal allelic effects and 
strict additivity, see Appendix S3. 



but see Rezende and Santos 2012). Therefore, it remains 
unclear how stochasticity might affect both the selection 
differential during selection for increased heat tolerance 
and/or the power to detect any eventual evolutionary 
response. Here, we employ computer simulations to show 
that selection responses and realized heritabilities depend 
on methodological context for knockdown temperature 
(i.e., the estimator), but are essentially context- indepen- 
dent for CT max (i.e., the parameter researchers attempt to 
estimate). As knockdown temperature involves a 



substantial amount of noise due to stochasticity, some 
ramping protocols may misleadingly suggest low evolu- 
tionary potential in CT max when genetic adaptation has in 
fact taken place. 

Computer Simulations 

We used computer simulations that mimic artificial selec- 
tion experiments for increasing knockdown temperature 
in the model species D. tnelanogaster, employing the 
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"Gompertz" code provided in Santos et al. (2011) to 
compute knockdown temperatures. The metric trait under 
selection was heat tolerance, measured as knockdown 
temperature in a ramping protocol. The target variable 
knockdown temperature presumably estimates the under- 
lying parameter CT max , as discussed in Santos et al. 
(2011), which is the polygenic character that we simulate. 
We assumed that short-term acclimation responses 
("hardening") did not occur during the thermal tolerance 
assays, which is an important caveat to keep in mind in 
experiments. 

In the initial model, CT max was determined by an 
arbitrary number of autosomal diallelic loci with purely 
additive effects for simplicity (model 1 below). Tolerance 
to high temperatures in large outbred populations of D. 
melanogaster is known to be a polygenic character, 
presumably affected by hundreds of genes on all chro- 
mosomes (Cavicchi et al. 1995; Loeschcke et al. 1997; 
Sorensen et al. 2005). However, one or a few candidate 
genes seem to explain much of the quantitative variation 
that we see in nature for heat knockdown temperature 
(Gilchrist and Huey 1999; Rand et al. 2010), which 
suggests that the distribution of allele effects affecting 
the quantitative variability in heat tolerance can have 
an exponential or geometric shape with increasingly 
fewer genes of progressively larger effects (Shrimpton 
and Robertson 1988; Orr 1999; Hayes and Goddard 
2001; but see Rockman 2012). As allele frequency 
changes during selection are highly dependent on allele 
frequencies and the distribution of allele effects, we also 
investigated the gamma distribution with shape parame- 
ter one-half and scale parameter one for the distribution 
of allele effects. 

Next, we incorporated additional genetic variation in 
metabolic rates (MR) to analyze how the selection proto- 
col might impact correlated responses on this trait that 
may, in turn, affect selection responses on heat tolerance 
(model 2). Water depletion during a thermal tolerance 
assay increases with metabolism (Rezende et al. 2011; see 
above), which is a serious concern in long assays that use 
slow heating rates because water content has a significant 
impact on heat tolerance (Maynard Smith 1957; Levins 
1969; Parsons 1980; Huey et al. 1992; Block et al. 1994). 
In our model, variation in MR affects the rate at which 
resources are depleted and, consequently, EC(ti) (see 
above). Recent debates emphasize the need to use ecologi- 
cally relevant slow heating rates when extrapolating 
laboratory estimates of heat tolerance to field conditions 
(Terblanche et al. 2007; Chown et al. 2009, 2010; Hoff- 
mann 2010; Mitchell and Hoffmann 2010; Sgro et al. 
2010). In this case, we predict a correlated response in 
decreasing MR when selecting flies for increased knockdown 
temperature. 



Model 1: Genetic variation for CT max 

The genotypic value for CT max had maximum range from 
a lower CT max (low) to an upper CT max (up) defined lim- 
its, fixed without loss of generality at 38°C and 42°C, 
respectively. We assumed CT max to be determined by 
/autosomal loci with two alleles each, 1 and 0 contrib- 
uting a = (42-38)/2/°C and 0°C, respectively; that is, we 
modeled purely additive effects with each locus contribut- 
ing equally to CT max (see Feder 1996). The genotypic 
value of each individual was thus 38°C plus a times the 
number of 1 alleles. Following Bulmer (1976), the total 
genotypic variance V G in this model, calculated from the 
distribution of CT max genotypic values, can be conve- 
niently partitioned into its three causal components: 

V G = V g + Dn- W + D L , (2) 

where V g is the genie variance computed from the 
observed allele frequencies; D H -w is the variance arising 
from deviations of perfect Hardy-Weinberg proportions 
at each locus; and D L is the variance due to deviations of 
linkage equilibrium among loci. Expressions for these 
variances are as follows (after Bulmer 1976): 

i 

V g = 2a 2 j2p.(l-Pi), 

i=i 

D H -w = a 1 [2Pj(0)Pi(2) - V2^(l)] , (3) 
i=i 

D L = 2a 2 ^cov(;\j), 

«<; 

where pi is the frequency of allele 1 at the rth locus; _?,(•) 
is the frequency of genotypes with 0, 1, and 2 alleles 1 at 
the z'th locus; and cov(ij) is the covariance between the 
number of alleles 1 at loci i and j (in the simulations D L 
is computed as D L = V G — V s — D H -w)- Assuming ran- 
dom mating D H -w varies around zero due to sampling 
fluctuations, but D L will depart from zero due to the 
effects of selection (Felsenstein 1965; Bulmer 1971, 1976). 
The / loci controlling CT max were assumed to be on the 
same chromosome for simplicity, and the recombination 
fraction between adjacent loci was zero in males (as it 
happens in Drosophila) and r in females with no interfer- 
ence. The recombination process followed the stochastic 
multilocus method described in Fraser and Burnell 
(1970). Briefly, the simulation of recombination involves 
the equivalent of a random walk along the length of the 
pair of homologous chromosomes, changing from one 
homologue to the other within the constraint of the 
probability of recombination. 

The phenotypic values for CT max were obtained by 
adding a normally distributed environmental component 
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with zero mean and variance V e to the genotypic values. 
Importantly, this environmental variance V e is the varia- 
tion in phenotype that cannot be explained by identifiable 
genetic differences and can be assumed to arise from 
uncontrolled random developmental variation among 
individuals. It affects the genetically determined basal 
CT max per se, and therefore has nothing to do with the 
additional residual variation that arises from stochasticity 
when estimating knockdown temperatures in experiments 
(Santos et al. 2011; see above). In other words, V e is the 
"real" environmental variance component that quantita- 
tive geneticists routinely introduce when modeling a 
metric trait (Falconer and Mackay 1996). 

The census population at each generation was N, with 
an equal number of females and males. All individuals 
were scored for their knockdown temperatures using the 
"Gompertz" computer code provided in Santos et al. 
(2011), and directional selection was applied by retaining 
the top 20% individuals from each sex at each generation 
(simulated individuals do not become sterilized after high 
temperature exposure; see Discussion). These selected 
individuals produced 2N gametes that were paired at ran- 
dom to render the next generation. The simulations were 
continued for g generations of directional selection. All 
simulations for model 1 assumed an average fruit fly 
weighing 1 mg with constant MR of 4.2 mL 0 2 /g/h at 18°C, 
or 0.07 /(L0 2 per fly min (Berrigan and Partridge 1997). 
Therefore, MR does not contribute to the residual varia- 
tion in CT max as assumed by Rezende et al. (2011). We 
also assumed that its total energy budget before the heat 
knockdown assay was equal to 171.6 jiL0 2 , and Q 10 = 3.5 
(see Santos et al. 2011). 

Model 2: Genetic variation for CT max and MR 

Our second model explicitly takes into account genetic 
differences in MR (known to be responsive to laboratory 
selection in Drosophila; Williams et al. 1997) to analyze 
how correlated responses in this trait could influence 
selection responses on heat tolerance. We recall that 
mortality rates in a heat tolerance assay may be partly 
determined by MR because this variable determines how 
fast water and energy resources are depleted (eq. 1; see 
also Rezende et al. 2011; Santos et al. 2011). We assumed 
/autosomal diallelic loci controlling MR with equal and 
additive effects for simplicity (see Simulation Results). 
The locations of all /+ /loci for CT max and MR were 
randomly assigned to the chromosome, and recombination 
frequency between adjacent loci was modeled as previ- 
ously indicated. 

The genotypic value for MR had maximum range from 
3.4 mL 0 2 /g/h to 5.6 mL 0 2 /g/h at 18°C (0.057 ^L0 2 
per fly min and 0.093 /(L0 2 per fly min, respectively). 



The reason for this range is that the initial frequencies of 
1 alleles increasing MR were randomly generated from a 
uniform distribution between 0.1 and 0.5, resulting in an 
average MR in the base population before selection close 
to 4.2 mL 0 2 /g/h at 18°C (0.07 /.(L0 2 per fly min) as in 
simulation model 1. The phenotypic values for MR were 
obtained by adding a normally distributed environmental 
component with zero mean and variance V' e to the geno- 
typic values. The average fruit fly also weighted lmg and 
its total energy budget prior to the heat knockdown assay 
was 171.6 ^L0 2 , and Q 10 = 3.5. 

The metric trait under selection was heat tolerance, 
measured as knockdown temperature as before. Note that 
any correlation arisen between CT max and MR during 
selection is not due to pleiotropy because the loci were 
assumed to affect each trait independently. Although link- 
age can be a cause of transient correlation, correlated 
responses in MR when up-selecting flies for knockdown 
temperature will be mainly caused by the "environment" 
(Falconer and Mackay 1996), which here means the 
methodology employed to score the flies (i.e., heating rate 
in the ramping assay). 

The simulation programs were implemented in the 
MATLAB algebra program environment (V7; Math Works 
2005) together with the collection of tools supplied by the 
Statistics Toolbox. The computer code is available upon 
request from the corresponding author. 

Simulation Results 

In all simulations, the census size was N = 5000 flies at 
each generation. Females and males were selected sepa- 
rately for knockdown temperature, and the top 20% of 
each sex was retained (i.e., 500 females and 500 males). 
The number of loci was /= 20, 40 with initial frequencies 
of alleles 1 randomly generated from a uniform distribu- 
tion between 0.1 and 0.5. Recombination frequencies were 
assumed to be r = 0.05, 0.15, 0.25 in females (population 
recombination frequency should be multiplied by one-half 
because Drosophila males lack recombination). Unless 
otherwise stated, the heritability of CT max (and MR; 
simulation model 2) was assumed to be h 2 = 0.25 in the 
base population before selection, a reasonable value for 
physiological traits (Mousseau and Roff 1987; Roff and 
Mousseau 1987). 

Model 1: Genetic variation for CT max 

Results from sample simulations with /= 20 and r = 0.25 
are plotted in Fig. 1A-D for flies selected with a fast 
ramping assay (T 0 = 28°C, AT = 0.5°C/min), and in 
Fig. 2A-D for flies selected with a slow ramping protocol 
(T 0 = 28°C, AT = 0.06°C/min). In both simulations, we 
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Figure 2. Sample numerical results from simulation model 1 using the same base population than in Fig. 1. Flies were subjected to 50 
generations of up-selection for knockdown temperature (the top 20% of each sex was retained) using a slow ramping protocol with To = 28°C 
and A7"= 0.06°C/min. Panel A plots the increase in knockdown temperature and its realized heritability, estimated by regressing the response to 
selection against the cumulated selection differential over the first 12 generations of selection. Panel B plots the increase in CT max , which was 2.5 
times higher than the increase in knockdown temperature. Panel C plots the frequency changes of alleles Y increasing CT ma> ., which eventually 
reached fixation. Panel D plots the total genotypic variance V G together with its causal components. Plots for CT max , allele frequencies and 
variance components are framed in shadow because their responses to selection for knockdown temperature are hidden to the experimentalist. 
For results of a similar model relaxing the assumptions of equal allelic effects and strict additivity, see Appendix S3. 



assumed additive and equal allele effects for CT max . The 
per-locus contribution to the additive genetic variance in 
the base population rang ed from 1.8 x 1(T 3 ( 0 C) 2 when 
p = 0.1 to 5.0 x 10" 3 (°C) 2 when p = 0.5. 

With fast ramping knockdown temperature rose from 
approximately 37°C to 39.5°C (Akt fast = 2.5°C), the real- 
ized heritability was 0.145 (Fig. 1A). However, under slow 
ramping, the increase in knockdown temperature was 
from approximately 36.2°C to 37.4°C (Akt slow = 1.2°C) 
with realized heritability 0.098 (Fig. 2A). In both situa- 
tions, the underlying heritability of CT max (h 2 « 0.25) 
was grossly underestimated. Importantly, the response of 
CT max to directional selection on knockdown temperature 
was essentially independent of the ramping conditions: it 
rose from an initial temperature of approximately 39°C 



and plateaued around the maximum attainable tempera- 
ture of 42°C (ACT max = 3°C) (Fig. IB, 2B) once alleles 
increasing CT max went to fixation (Fig. 1C, 2C). As 
expected from theory, the genotypic variance of CT max 
was reduced by directional selection because of the gener- 
ation of negative gametic linkage disequilibrium (Fig. ID, 
2D), the so-called "Buhner effect" (Bulmer 1971). 

We performed extensive computer simulations to cover 
a wide range of experimental conditions for all combina- 
tions of number of loci (f= 20, 40) determining CT max 
and recombination frequencies (r = 0.05, 0.15, 0.25). For 
each combination of parameter values, the same initial 
population of N = 5000 flies was subjected to 12 genera- 
tions of up-selection for knockdown resistance under 400 
different ramping protocols, after combining 20 initial 
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temperatures (T 0 ranging from 15°C to 34°C with interval 
1°C) with 20 heating rates (AT ranging from 0.05 to 1°C/ 
min with interval 0.05°C/min). Realized heritabilities for 
knockdown temperature, as well as the absolute increase 
of knockdown temperature and CT max after selection, 
against T 0 and AT are given as 3-D surface plots in 
Appendix S2. The numerical results provide a clear snap- 
shot on how different ramping rates affect the response of 
directional selection for knockdown temperature, the real- 
ized heritabilities, and the underlying genetic responses in 
CT max . Several conclusions emerge and can be summarized 
as follows: (1) The efficiency of selection in changing the 
allele frequencies and the trait mean (knockdown temper- 
ature, CT max ) obviously depends on the magnitude of 
allele effects, which decreases with an increasing number 
of loci £ (2) The effectiveness of selection is slightly lower 
when recombination frequencies are low (r = 0.05) as 
expected by the Hill-Robertson effect (Hill and Robertson 
1966; Felsenstein 1974), which establishes that the tighter 
is the linkage the higher is the perturbation that selection 
at one locus will have on other loci. (3) Realized herit- 
abilities for knockdown temperature are in the range 0.08 
-0.16 and grossly underestimate the true heritability of 
CT max (h 2 « 0.25), if both T 0 and AT are high realized 
heritabilities for knockdown temperature increase by 
about 60% when compared with slow ramping protocols 
that start with an initially low or moderate T 0 . (4) The 
increase of knockdown temperature(Akt) after 12 genera- 
tions of selection can differ up to two or threefold 
according to the methodology, with Akt being higher 
when both T 0 and AT are high. (5) The underlying 
increase in CT max is always higher than Akt and more or 
less independent of the methodological approach, with a 
maximum difference across ramping protocols generally 
less than 20%. 

Realized heritabilities for knockdown temperature will 
obviously change according to the underlying heritability 
of CT max . For instance, assuming h 2 = 0.1 for CT max in 
the base population simulations as those performed in 
Fig. 1A-D and Fig. 2A-D show that realized heritability 
for knockdown temperature is around 0.076 with fast 
ramping and 0.054 with slow ramping; and with h 2 = 0.4, 
the corresponding values are 0.186 with fast ramping and 
0.132 with slow ramping. However, the important point 
is that realized heritabilities for knockdown temperature 
will always underestimate the "true" heritability of CT max , 
and the underestimation is higher the slower the ramping 
rate. Needless to say, if the true heritability of CT max is 
low, the power to detect a realized heritability for knock- 
down temperature significantly different from zero under 
"ecologically realistic" slow ramping rates will likely be 
very low. What our simulations emphasize is that, with 
certain experimental approaches, it is impossible to 
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discriminate if low heritabilities reflect a biological phe- 
nomenon or a measurement artifact. 

Our next step was to analyze to what extent the previ- 
ous conclusions are robust to simplifying genetic details; 
namely, additive and equal allele effects for CT max . 
Assuming unequal allele effects with nonadditive contri- 
butions to CT max simulations show that the previous 
conclusions quantitatively hold (Appendix S3). 

Model 2: Incorporating genetic variation in 
metabolic rates 

As correlated responses in decreasing MR may occur 
when selecting flies for increased knockdown temperature, 
our final model explicitly takes into account genetic dif- 
ferences for both CT max and MR to analyze to what 
extent a reduction in MR can affect our previous conclu- 
sions, and also because variation in MR introduces an 
additional source of residual variance when scoring flies 
for knockdown temperature (see above). We assumed a 
simple additive and equal allele effects model because our 
previous simulations showed that numerical results were 
robust to these simplifying assumptions (Appendix S3). 

Results from sample simulations with /= 20 loci for 
both CT max and MR with r = 0.25 indicate that when 
flies are selected with a fast ramping assay (T 0 = 28°C, 
AT = 0.5°C/min), the loci affecting MR fluctuated more 
or less randomly and flies' average MR did not 
substantially change during selection (Fig. 3E, 3F). Inter- 
estingly, D. melanogaster flies selected for knockdown 
temperature under a fast heating rate (T 0 = 30°C, 
AT « 0.4°C/min) did not show correlated responses in 
MR in the upper thermal range (Table 1 in Folk et al. 
2007). Conversely, in the slow ramping protocol (T 0 = 
28°C, AT = 0.06°C/min), simulations show that there was 
a clear declining trend in the frequencies of alleles that 
raise MR (alleles 1 ): average MR at 18°C decreased from 
0.067 fj.L0 2 per fly min to 0.058 ^L0 2 per fly min after 
selection (Fig. 4E, 4F). 

However, these correlated responses did not have any 
major impact on how different ramping rates affect the 
realized heritabilities for knockdown temperature (0.152 
with fast ramping and 0.113 with slow ramping), 
although the increase in knockdown temperature after 
selection was approximately twice as higher with fast 
ramping (Akt fast = 2.5°C) than with slow ramping 
(Akt slow = 1.4°C) (Fig. 3A, 4A). The underlying genetic 
response in CT max was essentially the same in both cases: 
it rose up to 42 °C after 30 generations (Fig. 3B, 4B) once 
alleles increasing CT max went to fixation (Fig. 3C, 4C). As 
before, the genotypic variance in CT max was reduced by 
directional selection because of the generation of negative 
gametic linkage disequilibrium (Fig. 3D, 4D). Extensive 
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Figure 3. Sample numerical results from simulation model 2 assuming genetic variation for CT max and MR. The census size at each generation 
was W = 5, 000 flies. They were subjected to 50 generations of up-selection for knockdown temperature (the top 20% of each sex was retained) 
using a fast ramping protocol with To = 28°C and AT = 0.5°C/min. Both CT max and MR were independently controlled by 20 diallelic loci on the 
same chromosome, with recombination frequency between adjacent loci r= 0.25 in females. Allelic frequencies in the base population ranged 
from p = 0.1 to p = 0.5 for allele Y, which has additive effects 0.1°C for CT max or 0.055 mL0 2 /g/h at 18°C for MR. Heritabilities of CT max and 
MR were h 2 = 0.2B in the base population. Panel A plots the increase in knockdown temperature and its realized heritability, estimated by 
regressing the response to selection against the cumulated selection differential over the first 12 generations of selection. Panel B plots the 
increase in CT max , which was 20% higher relative to the increase in knockdown temperature. Panel C plots the frequency changes of alleles Y 
increasing CT ma>! , which eventually reached fixation. Panel D plots the total genotypic variance V G together with its causal components. The genie 
variance component V g initially increased (approximately up to generation 8) as a consequence of the changes in allele frequencies, but the 
genotypic variance V G steadily decreased because of linkage disequilibrium (D L ). Panel E plots the correlated response in MR, which did not 
substantially change during selection. Panel F plots the frequencies of alleles Y increasing MR, which fluctuated more or less randomly during 
selection. Plots for CT max , allele frequencies, and variance components are framed in shadow because their responses to selection for knockdown 
temperature are hidden to the experimentalist. 
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Figure 4. Sample numerical results from simulation model 2 using the same base population than in Fig. 3. Flies were subjected to 50 
generations of up-selection for knockdown temperature (the top 20% of each sex was retained) using a slow ramping protocol with T 0 = 28°C 
and A7"= 0.06°C/min. Panel A plots the realized heritability of knockdown temperature estimated by regressing the response to selection against 
the cumulated selection differential over the first 12 generations of selection. Panel B plots the increase in CT max alter selection, which was 115% 
higher relative to the increase in knockdown temperature. Panel C plots the frequency changes of alleles T increasing CT max , which eventually 
reached fixation. Panel D plots the total genotypic variance V G together with its causal components. Panel E plots the correlated response in MR, 
which dropped by approximately 16% after selection for increasing knockdown temperature. Panel F plots the frequency changes of alleles Y 
affecting MR, which clearly decreased during selection and were eventually lost. Plots for CT max , variance components, and allele frequencies are 
framed in shadow because their responses to selection for knockdown temperature are hidden to the experimentalist. 
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computer simulations as those performed in Appendix S2 
reinforce these conclusions (results not shown). 

Finally, some simulated scenarios (Appendix S2) can be 
criticized for being highly unrealistic in practical terms. If 
the initial temperature in the ramping assay is low, say 
To < 20°C, the time taken to score the flies for knock- 
down temperature at any generation can be higher than 
3 h with slow heating rates, and acclimation effects dur- 
ing the assays cannot probably be ignored (Sgro et al. 
2010; Santos et al. 2011). However, for any initial temper- 
ature researchers might consider appropriate to perform 
the experiment (e.g., T 0 = 28°C as used in the sample 
simulations), the previous conclusions hold. 

Discussion 

A number of authors have recently questioned the stan- 
dard belief that there is abundant genetic variation in any 
ecologically relevant trait for natural selection to act on 
(Blows and Hoffmann 2005). Some tropical rainforest 
Drosophila species appear to have lost heritable variation 
for desiccation resistance (Hoffmann et al. 2003a; Keller- 
mann et al. 2006, 2009), an important physiological trait 
that might have an impact on species distributions. Selection 
experiments for increasing heat knockdown resistance in 
the cosmopolitan species D. melanogaster also suggest low 
but significant levels of genetic variation, with realized 
heritabilities around 7-12% (McColl et al. 1996; Bubli 
et al. 1998; Gilchrist and Huey 1999). However, flies in 
the experiments were acutely exposed to a high tempera- 
ture or to a fast heating rate, which has raised concerns 
about extrapolations to field conditions because estimates 
under "ecologically realistic" slow heating rates suggest 
that the heritability of upper thermal limits is close to 
zero (Mitchell and Hoffmann 2010). These findings add 
considerable fuel to the debate on whether species will be 
able to persist by adapting genetically to our changing cli- 
mate (Parmesan 2006; Skelly et al. 2007; Visser 2008; 
Chevin et al. 2010). 

The speculation that the heritability of CT max can 
change as a direct consequence of the heating rate was 
advanced by Chown et al. (2009). For instance, Mitchell 
and Hoffmann (2010) showed that the estimated herita- 
bility of heat tolerance in one Australian population 
(Gordonvale) of D. melanogaster dropped from 
h 2 = 0.22 ± 0.07 when flies were placed acutely at a 
stressful temperature of 38°C to h 2 = 0.05 ± 0.07 when 
assayed with a slow heating rate of 0.06°C/min. Both 
Chown et al. (2009) and Mitchell and Hoffmann (2010) 
explicitly talked about CT max , although what they mea- 
sured was heat knockdown resistance using different 
methods. We have cautioned on the common misconcep- 
tion of equating CT max with knockdown resistance 



(Santos et al. 2011; see above), and the present results 
dramatically illustrate the repercussions of our warnings. 
The value h 2 = 0.05 ± 0.07 in Mitchell and Hoffmann 
(2010) is a heritability estimate of knockdown tempera- 
ture under slow ramping. It is entirely consistent to 
simultaneously have a moderate heritability of CT max and 
a low heritability of knockdown temperature in dynamic 
experiments that use a ramping rate of 0.06°C/min 
(Appendix S2). Their conclusion that upper thermal 
limits have low evolutionary potential under ecologically 
relevant slow heating rates is consequently incorrect. The 
estimate h 2 = 0.22 ± 0.07 (Table 2 in Mitchell and Hoff- 
mann 2010) is probably closer to the true heritability of 
CT max in their population (qualitatively similar results 
were reported for a Melbourne population; see Mitchell 
and Hoffmann 2010). 

Why selection responses are context- 
dependent for knockdown temperature, but 
not for CT max ? 

The amount of noise introduced by stochasticity during a 
heat tolerance assay obviously lowers the correlation 
between knockdown temperature (estimator) and CT max 
(parameter), and we have shown that the expected repeat- 
ability of heat tolerance can be very low (Santos et al. 
2011). Actually, the limited evidence available in Drosoph- 
ila points to a repeatability of 20% (Krebs and Loeschcke 
1997; see also Santos et al. 2011). However, when select- 
ing for increasing heat tolerance, the crucial point is to 
know how the selection differential applied to this pheno- 
typic character translates to the genetic differences in the 
causal trait CT max . In our theoretical approach to model 
heat resistance assays (eq. 1), we explicitly incorporated 
the obvious assumption that survival probability steadily 
decreases toward zero when body temperature T b -> 
CT max . Therefore, it may happen that an individual that 
can tolerate, in the very best of cases, a temperature of 
41°C collapses at a lower T b of, say, 39°C because the 
time-dependent cumulative probability of dying at this 
temperature is higher than zero (importantly, the lower 
the heating rate in a ramping assay, the higher the proba- 
bility of dying at a temperature substantially lower than 
CT max ; Santos et al. 2011). But, this cannot happen the 
other way around. This asymmetry guarantees that the 
top percentile of selected individuals (20% in the simula- 
tions) for knockdown temperature includes most of the 
"best" individuals for CT max , and this is largely independent 
of heating rates. For instance, assume a computer-generated 
base population from simulation model 1 selected for 
knockdown temperature. A selection differential of 0.91°C 
(top 20%) under fast ramping translates to a selection dif- 
ferential of 0.67°C for CT max , and a selection differential 
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of 0.39°C (top 20%) under slow ramping to a selection 
differential of 0.63°C for CT max . Although the exact mag- 
nitude of these numbers can obviously change according 
to the specific details of the function that describes the 
survival probability during a heat resistance assay, it is 
undoubtedly true that the previous asymmetry holds. 

Selection for increased heat tolerance 

Our results suggest that CT max can evolve substantially 
more in response to selection than is estimated by empiri- 
cal measures of knockdown temperature and, most 
importantly, this response seems to be independent of 
heating rates (correlated responses in MR, on the other 
hand, are more pronounced under slow heating rates; cf. 
Figs. 3E, 4E). Although computer simulations have 
focused in D. melanogaster for obvious reasons, the prob- 
lem is absolutely general because the estimation of upper 
thermal limits involves placing the individuals under 
stressful conditions and record the time to collapse or 
death. Thus, even though resource depletion during a 
heat resistance assay may not be a major concern for a 
larger organism than a Drosophila, stochasticity is an 
unavoidable source of error that downwardly bias esti- 
mates of heat tolerance and its evolutionary potential. 
This can be easily illustrated by setting EC(t;) w 0 in our 
model (eq. 1); that is, by assuming that in the limit, the 
physiological condition of individuals does not decay during 
the heat tolerance assay. Simulations show that realized her- 
itabilities for knockdown temperature underestimate the 
"true" heritability of CT max , and the underestimation is 
higher the slower the ramping rate (results not shown). 

These findings have important repercussions for our 
understanding of the evolution of thermal tolerance for 
two reasons. First, they suggest that the methodology 
employed can seriously underestimate the evolutionary 
response of CT max in selection experiments, which is in 
close agreement with our findings of shallower latitudinal 
clines due to methodology (Santos et al. 2011). Selection 
with the "knockdown tube" (Huey et al. 1992) - an appa- 
ratus in which the temperature or time at which flies lose 
ability to cling on the walls of the tube and fall down can 
be readily recorded - has provided empirical evidence of 
the evolutionary potential of heat tolerance, resulting in 
an increase of nearly 2.5°C in knockdown temperature 
and a realized heritability of roughly 0.12 for this trait 
(Gilchrist and Huey 1999; see also McColl et al. 1996). 
Our results are in close agreement with these values, but 
suggest that the overall response to selection and 
heritability for the parameter CT max may be even higher 
(although simulations are not meant to mimic the 
conditions of the knockdown tube, a decrease in average 
heat tolerance is expected simply due to stochasticity). 



This theoretical framework can also explain why flies 
successfully selected for increasing knockdown resistance 
in the knockdown tube do not show a higher knockdown 
resistance than their respective controls when assayed in 
glass vials (Hoffmann et al. 1997). Although this discrep- 
ancy suggested to some researchers that the physiological 
and genetic mechanisms accounting for heat tolerance 
vary according to the methodology employed (Hoffmann 
et al. 1997, 2003b; Rako et al. 2007; Sgro et al. 2010), we 
showed that the absence of correlation between heat tolerance 
indices is not evidence of different underlying mechanisms 
(Santos et al. 2011). Nonetheless, gender-specific patterns 
show that indices of physiological tolerance differ between 
methods: whereas, in the knockdown tube, males are 
somewhat more resistant to knockdown than females (Hoff- 
mann et al. 1997; p. 394; see also Jenkins and Hoffmann 
1994; Bubli et al. 1998), a common result in glass vials is 
that D. melanogaster females are more resistant to heat 
stress than males (Huey et al. 1991; Loeschcke et al. 1997; 
Mitchell and Hoffmann 2010; Parkash et al. 2010). These 
methodological differences (see also Folk et al. 2006) 
should be critically addressed before speculating about a 
putative independent genetic control of heat tolerance 
indices inferred from correlated responses (or lack 
thereof) across methodologies. 

Second, our simulations show that the effects of sto- 
chasticity on mortality as temperatures increase do not 
have a major impact on the selection differential of 
CT max , hence genes for increased heat tolerance should be 
eventually selected regardless of the thermal regime. Thus, 
the main issue from an empirical perspective remains 
detecting, rather than eliciting, an evolutionary response 
(see also Santos et al. 2011). The prediction stemming 
from our results is that heat tolerance will increase to 
roughly the same level regardless the ramping rates 
employed during selection, which can be tested with the 
knockout tube or undertaking family selection experi- 
ments using the knockdown vial technique (which is 
advantageous because family means would be less affected 
by methodological noise, and selected flies would not be 
stressed). Similarly, should our results be extrapolated to 
natural conditions, they would suggest that daily and sea- 
sonal variations in the rate of change in temperature have 
only a minor effect on the overall evolutionary response 
in CT max , everything else being equal. 

Concluding remarks 

This study arose from the paradox that the most wide- 
spread and common Drosophila species (Powell 1997) 
apparently exhibits limited adaptive potential for upper 
thermal limits (Mitchell and Hoffmann 2010). This con- 
clusion is in striking conflict with the invasive success of 
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D. melanogaster (which originated in sub-Saharan Africa 
and established in Europe and Asia, more recently, in 
both the New World and Australia), a species subjected 
to spatially varying selection for many traits, including 
thermotolerance (e.g., David and Capy 1988; Sezgin et al. 
2004; Schmidt et al. 2005; Hoffmann and Weeks 2007; 
Gonzalez et al. 2010). Here, we demonstrate that this 
contradiction may stem from the confusion between 
parameter (CT max ) and estimator (knockdown temperature), 
with the unfortunate result that "ecologically realistic" 
assays yield highly downwardly biased estimates of upper 
thermal limits (Rezende et al. 2011; Santos et al. 2011) 
and their true evolutionary potential. A reviewer of San- 
tos et al.'s (2011) paper considered the take-home mes- 
sage that estimates of CT max are highly sensitive to 
methodology as rather depressing. Ironically, this message 
turns out to be good news because adaptive genetic 
responses for increasing upper thermal limits may be 
higher than acknowledged in recent studies. 

Our results also illustrate how the experimental 
approaches adopted might substantially affect the conclu- 
sions drawn from a particular investigation (Chown et al. 
2009; p. 138). Importantly, the solution to this problem 
does not entail compiling heat tolerance estimates under 
a myriad of conditions and increasingly intricate experi- 
ments (e.g., Terblanche et al. 2008; Chidawanyika and 
Terblanche 2011; Overgaard et al. 2011a). It is currently 
clear that the uncontrolled effects of cumulative thermal 
stress, its impact on correlated traits and on intrinsic sur- 
vival probabilities may have a larger impact on heat toler- 
ance estimates than the factors under study (Rezende 
et al. 2011; Santos et al. 2011). Many of the indices 
employed in recent studies are hardly comparable across 
systems, their precision and validity have not been 
assessed experimentally (repeatability estimates of mea- 
surements of thermal tolerance and other physiological 
limits are virtually absent; Santos et al. 2011; Krebs and 
Loeschcke 1997; see also Wolak et al. 2011; Castaneda 
et al. 2012) and neither has their "ecological relevance" 
(Rezende and Santos 2012). Therefore, resulting patterns 
should be carefully assessed in the light of our findings, 
which clearly show that there is substantially more to 
thermal limits and their potential to respond to selection 
than meets the eye. 

Acknowledgments 

We thank five anonymous reviewers for their insightful 
comments on earlier drafts that helped to significantly 
improve the manuscript. M.S. is supported by grant 
CGL2010-15395 from the Ministerio de Ciencia 
e Innovacion (Spain) and by the ICREA Academia program. 
L.E.C. is supported by a Juan de la Cierva fellowship 



(JCI-2010-06156) from the Ministerio de Ciencia e Inno- 
vacion. E.L.R. is supported by a Ramon y Cajal contract 
and by grant BFU2009-07564 from the Ministerio de 
Ciencia e Innovacion. Financial support by grant 
2009SGR 636 from Generalitat de Catalunya to the Grup 
de Biologia Evolutiva is also gratefully acknowledged. 

Conflict of Interest 

None declared. 
References 

Addo-Bediako, A., S. L. Chown, and K. J. Gaston. 2000. 

Thermal tolerance, climatic variability and latitude. Proc. 

R. Soc. London B 267:739-745. 
Angilletta, M. J., Jr. 2009. Thermal adaptation. A theoretical 

and empirical synthesis. Oxford Univ. Press, Oxford. 
Berrigan, D., and L. Partridge. 1997. Influence of temperature 

and activity on the metabolic rate of adult Drosophila 

melanogaster. Comp. Bioch. Physiol. 118A:1301-1307. 
Block, W., N. R. Webb, S. Coulson, I. D. Hodkinson, and 

M. R. Worland. 1994. Thermal adaptation in the Arctic 

collembolan Onychinrus arcticus (Tullberg). J. Insect Physiol. 

40:715-722. 

Blows, M. W., and A. A. Hoffmann. 2005. A reassessment of 
genetic limits to evolutionary change. Ecology 86:1371-1384. 

Bubli, O. A., A. G. Imasheva, and V. Loeschcke. 1998. 
Selection for knockdown resistance to heat in Drosophila 
melanogaster at high and low larval densities. Evolution 
52:619-625. 

Bulmer, M. G. 1971. The effect of selection on genetic 

variability. Am. Nat. 105:201-211. 
Bulmer, M. G. 1976. The effect of selection on genetic 

variability: a simulation study. Genet. Res. 28:101-117. 
Castaneda, L. E., G. Calabria, L. A. Betancourt, E. L. Rezende, 

and M. Santos. 2012. Measurement error in heat tolerance 

assays. J. Therm. Biol 37:432-437. 
Cavicchi, S., D. Guerra, V. Latorre, and R. B. Huey. 1995. 

Chromosomal analysis of heat-shock tolerance in Drosophila 

melanogaster evolving at different temperatures in the 

laboratory. Evolution 49:676-684. 
Chevin, L.-M., R. Lande, and G. M. Mace. 2010. Adaptation, 

plasticity, and extinction in a changing environment: 

towards a predictive theory. PLoS Biol. 8(4):el000357. 
Chidawanyika, F., and J. S. Terblanche. 2011. Rapid thermal 

responses and thermal tolerance in adult codling moth 

Cydia pomonella (Lepidoptera: Tortricidae). J. Insect 

Physiol. 57:108-117. 
Chown, S. L., K. R. Jumbam, J. G. Sorensen, and J. S. 

Terblanche. 2009. Phenotypic variance, plasticity 

and heritability estimates of critical thermal limits 

depend on methodological context. Funct. Ecol. 23: 

133-140. 



2012 The Authors. Published by Blackwell Publishing Ltd. 



2877 



Evolution of Upper Thermal Limits 



Santos et al. 



Chown, S. L., A. A. Hoffmann, T. N. Kristensen, M. J. 

Angilletta, N. C. Stenseth, and C. Pertoldi. 2010. Adapting 
to climate change: a perspective from evolutionary 
physiology. Clim. Res. 43:3-15. 

David, J. R., and P. Capy. 1988. Genetic variation of Drosophila 
melanogaster natural populations. Trends Genet. 4:106-111. 

Deutsch, C. A., J. J. Tewksbury, R. B. Huey, K. S. Sheldon, 
C. K. Ghalambor, D. C. Haak, et al. 2008. Impacts of 
climate warming on terrestrial ectotherms across latitude. 
Proc. Natl. Acad. Sci. USA 105:6668-6672. 

Elliott, J. M., J. A. Elliott, and J. D. Allonby. 1994. The critical 
thermal limits for the stone loach, Noemacheilus barbatulus, 
from three populations in north-west England. Freshwater 
Biol. 32:593-601. 

Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to 
quantitative genetics. Longman Group Limited, Harlow. 

Feder, M. E. 1996. Ecological and evolutionary physiology of 
stress proteins and the stress response: the Drosophila 
melanogaster model. Pp. 79-102 in I. A. Johnston and 
A. F. Bennet, eds. Phenotypic and evolutionary adaptation 
to temperature. Cambridge Univ. Press, Cambridge. 

Felsenstein, J. 1965. The effect of linkage on directional 
selection. Genetics 52:349-363. 

Felsenstein, J. 1974. The evolutionary advantage of 
recombination. Genetics 78:737-756. 

Folk, D. G., P. Zwollo, D. M. Rand, and G. W. Gilchrist. 2006. 
Selection on knockdown performance in Drosophila 
melanogaster impacts thermotolerance and heat-shock 
response differently in females and males. J. Exp. Biol. 
209:3964-3973. 

Folk, D. G., L. A. Hoekstra, and G. W. Gilchrist. 2007. Critical 

thermal maxima in knockdown-selected Drosophila: are 

thermal endpoints correlated? J. Exp. Biol. 210:2649-2656. 
Fraser, A., and D. Burnell. 1970. Computer models in genetics. 

McGraw-Hill, New York. 
Ghalambor, C. K., J. K. McKay, S. P. Carroll, and D. N. 

Reznick. 2007. Adaptive versus non-adaptive phenotypic 

plasticity and the potential for contemporary adaptation in 

new environments. Funct. Ecol. 21:394-407. 
Gilchrist, G. W., and R. B. Huey. 1999. The direct response of 

Drosophila melanogaster to selection on knockdown 

temperature. Heredity 83:15-29. 
Gonzalez, J., T. L. Karasov, P. W. Messer, and D. A. Petrov. 

2010. Genome-wide patterns of adaptation to temperate 

environments associated with transposable elements in 

Drosophila. PLoS Genet. 6(4):el000905. 
Hayes, B., and M. E. Goddard. 2001. The distribution of the 

effects of genes affecting quantitative traits in livestock. 

Genet. Sel. Evol. 33:209-229. 
Helmuth, B., J. G. Kingsolver, and E. Carrington. 2005. 

Biophysics, physiological ecology, and climate change: does 

mechanism matter? Ann. Rev. Physiol. 67:177-201. 
Hill, W. G., and A. Robertson. 1966. The effect of linkage on 

limits to artificial selection. Genet. Res. 8:269-294. 



Hochachka, P. W., and G. N. Somero. 2002. Biochemical 

adaptation: mechanism and process in physiological 

evolution. Oxford Univ. Press, New York. 
Hoffmann, A. A. 2010. Physiological climatic limits in 

Drosophila: patterns and implications. J. Exp. Biol. 

213:870-880. 

Hoffmann, A. A., and C. M. Sgro. 2011. Climate change and 
evolutionary adaptation. Nature 470:479-485. 

Hoffmann, A. A., and A. R. Weeks. 2007. Climatic selection 
on genes and traits after a 100 year-old invasion: a critical 
look at the temperate-tropical clines in Drosophila 
melanogaster from eastern Australia. Genetica 129:133-147. 

Hoffmann, A. A., H. Dagher, M. Hercus, and D. Berrigan. 
1997. Comparing different measures of heat resistance in 
selected lines of Drosophila melanogaster. J. Insect Physiol. 
43:393-405. 

Hoffmann, A. A., R. J. Hallas, J. A. Dean, and M. Schiffer. 
2003a. Low potential for climatic stress adaptation in a 
rainforest Drosophila species. Science 301:100-102. 

Hoffmann, A. A., J. G. S0rensen, and V. Loeschcke. 2003b. 
Adaptation of Drosophila to temperature extremes: bringing 
together quantitative and molecular approaches. J. Therm. 
Biol 28:175-216. 

Huey, R. B., L. Partridge, and K. Fowler. 1991. Thermal 

sensitivity of Drosophila melanogaster responds to laboratory 
natural selection. Evolution 45:751-756. 

Huey, R. B., W. D. Crill, J. G. Kingsolver, and K. E. Weber. 
1992. A method for rapid measurement of heat or cold 
resistance of small insects. Funct. Ecol. 6:489-494. 

Huey, R. B., C. A. Deutsch, J. J. Tewksbury, L. J. Vitt, P. E. 
Hertz, H. J. Alvarez Perez, et al. 2009. Why tropical forest 
lizards are vulnerable to climate warming. Proc. R. Soc. 
London B 276:1939-1948. 

Jenkins, N. L., and A. A. Hoffmann. 1994. Genetic and 

maternal variation for heat resistance in Drosophila from the 
field. Genetics 137:783-789. 

Kellermann, V. M., B. van Heerwaarden, A. A. Hoffmann, and 
C. M. Sgro. 2006. Very low additive genetic variance and 
evolutionary potential in multiple populations of two 
rainforest Drosophila species. Evolution 60:1104-1108. 

Kellermann, V., B. van Heerwaarden, C. M. Sgro, and 

A. A. Hoffmann. 2009. Fundamental evolutionary limits in 
ecological traits drive Drosophila species distributions. 
Science 325:1244-1246. 

Krebs, R. A., and V. Loeschcke. 1997. Estimating heritability in 
a threshold trait: heat-shock tolerance in Drosophila buzzatii. 
Heredity 79:252-259. 

Levins, R. 1969. Thermal acclimation and heat resistance in 
Drosophila species. Am. Nat. 103:483-499. 

Loeschcke, V., R. A. Krebs, J. Dahlgaard, and P. Michalak. 
1997. High temperature stress and the evolution of thermal 
resistance in Drosophila. Pp. 175-190 in R. Bijlsma and 
V. Loeschcke, eds. Environmental stress adaptation and 
evolution. Birkhauser Verlag, Basel. 



2878 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



Santos et al. 



Evolution of Upper Thermal Limits 



Lutterschmidt, W. I., and V. H. Hutchinson. 1997. The critical 
thermal maximum: data to support the onset of spasms as 
the definitive end point. Can. J. Zool. 75:1553-1560. 

MathWorks. 2005. MATLAB. Version 7.0.4. MathWorks, Inc., 
Natick, MA. Available via <http://www.mathworks.com>. 

Maynard Smith, J. 1957. Temperature tolerance and 

acclimatization in Drosophila subobscura. J. Exp. Biol. 34:85-96. 

McColl, G., A. A. Hoffmann, and S. W. McKechnie. 1996. 
Response of two heat shock genes to selection for 
knockdown heat resistance in Drosophila melanogaster. 
Genetics 143:1615-1627. 

Mitchell, K. A., and A. A. Hoffmann. 2010. Thermal ramping rate 
influences evolutionary potential and species differences for 
upper thermal limits in Drosophila. Funct. Ecol. 24:694-700. 

Mora, C., and M. F. Maya. 2006. Effect of the rate of 
temperature increase of the dynamic method on the heat 
tolerance of fishes. J. Therm. Biol 31:337-341. 

Mousseau, T. A., and D. A. Roff. 1987. Natural selection and 
the heritability of fitness components. Heredity 59:181-197. 

Orr, H. A. 1999. The evolutionary genetics of adaptation: a 
simulation study. Genet. Res. 74:207-214. 

Overgaard, J., A. A. Hoffmann, and T. N. Kristensen. 2011a. 
Assesing population and environmental effects on thermal 
resistance in Drosophila melanogaster using ecologically 
relevant assays. J. Therm. Biol 36:409-416. 

Overgaard, J., T. N. Kristensen, K. A. Mitchell, and 

A. A. Hoffmann. 2011b. Thermal tolerance and tropical 
Drosophila species: does phenotypic plasticity increase with 
latitude? Am. Nat. 178:S80-S96. 

Parkash, P., V. Sharma, and B. Kalra. 2010. Correlated 
changes in thermotolerance traits and body color 
phenotypes in montane populations of Drosophila 
melanogaster: analysis of within- and between-population 
variations. J. Zool. 280:49-59. 

Parmesan, C. 2006. Ecological and evolutionary responses to 
recent climate change. Annu. Rev. Ecol. Evol. Syst. 37:637-669. 

Parsons, P. A. 1980. Parallel climatic races for tolerances to 
high temperature-desiccation stress in two Drosophila 
species. J. Biogeog. 7:97-101. 

Peck, L. S., M. S. Clark, S. A. Morley, A. Massey, and 

H. Rossetti. 2009. Animal temperature limits and ecological 
relevance: effects of size, activity and rates of change. Funct. 
Ecol. 23:248-256. 

Powell, J. R. 1997. Progress and prospects in evolutionary biology. 
The Drosophila model. Oxford Univ. Press, New York. 

Rako, L., M. J. Blacket, S. W. McKechnie, and A. A. 

Hoffmann. 2007. Candidate genes and thermal phenotypes: 
identifying ecologically important genetic variation for 
thermotolerance in the Australian Drosophila melanogaster 
cline. Mol. Ecol. 16:2948-2957. 

Rand, D. M., D. M. Weinreich, D. Lerman, D. Folk, and G. 
W. Gilchrist. 2010. Three selections are better than one: 



clinal variation of thermal QTL from independent 
selection experiments in Drosophila. Evolution 64:2921- 
2934. 

Rezende, E. L., and M. Santos. 2012. Comment on 

'Ecologically relevant measures of tolerance to potentially 
lethal temperatures'. J. Exp. Biol. 215:702-703. 

Rezende, E. L., M. Tejedo, and M. Santos. 2011. Estimating 
the adaptive potential of critical thermal limits: 
methodological problems and evolutionary implications. 
Funct. Ecol. 25:111-121. 

Ribeiro, P. L., A. Camacho, and C. A. Navas. 2012. 

Considerations for assessing maximum critical temperatures 
in small ectothermic animals: insights from leaf-cutting ants. 
PLoS ONE 7:e32083. 

Rockman, M. V. 2012. The QTN program and the alleles that 
matter for evolution: all that's gold does not glitter. 
Evolution 66:1-17. 

Roff, D. A., and T. A. Mousseau. 1987. Quantitative genetics 
and fitness: lessons from Drosophila. Heredity 58:103-118. 

Santos, M., L. E. Castaneda, and E. L. Rezende. 2011. Making 
sense of heat tolerance estimates in ectotherms: lessons from 
Drosophila. Funct. Ecol. 25:1169-1180. 

Schmidt, P. S., L. Matzkin, M. L. Ippolito, and W. F. Eanes. 
2005. Geographic variation in diapause incidence, life 
history traits, and climatic adaptation in Drosophila 
melanogaster. Evolution 59:1721-1732. 

Sezgin, E., D. D. Duvernell, L. M. Matzkin, Y. Duan, C.-T. 
Zhu, B. C. Verrelli, et al. 2004. Single locus latitudinal clines 
in metabolic genes, derived alleles, and their relationship to 
temperate adaptation in Drosophila melanogaster. Genetics 
168:923-931. 

Sgro, C. M., J. Overgaard, T. N. Kristensen, K. A. Mitchell, F. E. 
Cockerell, and A. A. Hoffmann. 2010. A comprehensive 
assessment of geographic variation in heat tolerance and 
hardening capacity in populations of Drosophila melanogaster 
from eastern Australia. J. Evol. Biol. 23:2484-2493. 

Shrimpton, A. E., and A. Robertson. 1988. The isolation of 
polygenic factors controlling bristle score in Drosophila 
melanogaster. II. Distribution of third chromosome 
bristle effects within chromosome sections. Genetics 
118:445-459. 

Skelly, D. K, L. N. Joseph, H. P. Possingham, L. K. 

Freidenburg, T. J. Farrugia, M. T. Kinnison, et al. 2007. 
Evolutionary responses to climate change. Cons. Biol. 
21:1353-1355. 

Sorensen, J. G., M. M. Nielsen, M. Kruhoffer, J. Justesen, and 
V. Loeschcke. 2005. Full genome gene expression analysis of 
the heat stress response in Drosophila melanogaster. Cell 
Stress Chap. 10:312-328. 

Sunday, J. M., A. E. Bates, and N. K. Dulvy. 2011. Global 
analysis of thermal tolerance and latitude in ectotherms. 
Proc. R. Soc. London B 278:1823-1830. 



2012 The Authors. Published by Blackwell Publishing Ltd. 



2879 



Evolution of Upper Thermal Limits 



Santos et al. 



Terblanche, J. S., J. A. Deere, S. Clusella-Trullas, C. Janion, and S. 
L. Chown. 2007. Critical thermal limits depend on 
methodological context. Proc. R. Soc. London B 274:2935-2942. 

Terblanche, J. S., S. Clusella-Trullas, J. A. Deere, and S. L. 
Chown. 2008. Thermal tolerance in a south-east African 
population of the tsetse fly Glossina pallidipes (Diptera, 
Glossinidae): implications for forecasting climate change 
impacts. J. Insect Physiol. 54:114-127. 

Visser, M. E. 2008. Keeping up with a warming world; 
assessing the rate of adaptation to climate change. Proc. 
R. Soc. London B 275:649-659. 

Williams, A. E., M. R. Rose, and T. J. Bradley. 1997. CO z 
release patterns in Drosophila melanogaster. the effect of 
selection for desiccation resistance. J. Exp. Biol. 200:615-624. 

Wolak, M. E., D. J. Fairban, and Y. R. Paulsen. 2011. 

Guidelines for estimating repeatability. Methods Ecol. Evol. 
3:129-137. 



Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Appendix SI. Stochasticity effects in heat knockdown 
resistance assays. 

Appendix S2. Numerical results from simulation model 1 
assuming additive and equal allele effects for CT max . 
Appendix S3. Genetic Variation for CT max assuming 
unequal allele effects and nonadditivity. 

Please note: Wiley- Blackwell are not responsible for the 
content or functionality of any supporting materials sup- 
plied by the authors. Any queries (other than missing 
material) should be directed to the corresponding author 
for the article. 



2880 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



